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Summary: The bi-continuum model composed of two interpenetrating and dynamically coupled 
material continua is analysed as a simplified but relatively accurate way to describe some physical phenomena 
in crystalline solids. The essential novelty of our approach consists in treating a crystalline medium as a 
bi-continuum, even if the crystalline lattice is structurally single-component. Particular attention is paid to 
the oscillatory behaviour of solutions on the atomic level. Starting from a discrete atomic chain, the basic 
formulation of the bi-continuum model is derived. The essential features of the model, including accuracy of 
the results as functions of physical parameters, are discussed. 
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1 Introduction 



Discrete structural models such as atomic lattices or stratified structures are capable of 
delivering, as a rule, highly relialable and accurate results concerning a wide range of physical 
phenomena in solid matter. Most frequently, however, to produce such results, an atomic 
model has to be treated with the aid of computer simulation techniques, what imposes 
severe restrictions on the size of computed samples and leads to the lack of generality in 
the results. Moreover, atomic models are sensitive to details of the interatomic potentials; 
the final reliability of results depends critically on the reliability of the potentials used. In 
consequence, simplified models which may produce analytical - although less accurate - 
results are still highly desirable. The loss of accuracy should not, however, be catastrophic. 

In the present paper we focus our attention on atomic configurations in complex layered 
structures or in the vicinity of planar or quasi-planar crystalline interfaces M. In such 
situations, the configurations frequently exhibit a sign-alternating behaviour |], [| on atomic 
scale, not reproducible by a simple classic continuum. We have found, however, that such 
a behaviour can be described by an appropriately chosen two-component continuum model, 
called the bi-continuum. 

The concept of multi-continuum, understood as a material medium on which several 
displacement fields are defined, is intuitively evident when one considers multicomponent 
media (like e.g. mixture of substances with different physical properties or the medium 
composed of electrons considered as one component and heavy ions as the other one). 

The idea given in || was to introduce the multicontinuum description for structurally 
complex media, even if such a medium has only one component. The different displacement 
fields are defined on different substructures. We shall show how the idea of multicontinuum 
works in the case when multi=2. 

From the point of view of independent variables, the problem is essentially one-dimensional. 
The corresponding variable will be denoted by x with, according to the needs, appropriate 
indices or diacritics. For the sake of simplicity, we perform our theoretical constructions on 
linear atomic chains, although in the real three-dimensional applications the "atoms" are 
intended to represent whole crystalline planes or appropriate structural layers. 

In Sections 2 and 3 we state the problem for a homogeneous atomic chain and examine 
the sign-alternating solutions near the "interface". In Section 4 we analyse a composite 
linear chain, consisting of two components. This provides the basis for construction of the 
bi-continuum model in Section 5. Subsequently we solve the interface problem in atomic and 
bi-continuum models, respectively, compare the solutions, and discuss the resulting degree 
of approximation. 

2 Homogeneous atomic chain 

Let us consider an atomic chain with interactions between the first and second nearest neigh- 
bours. Let Ui denote the displacement of the atom number i from its equilibrium position. 
By Kf,n,f and k s , n s we denote the elastic constant and initial stretch connected with shorter 
and longer bonds, respectively. 
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Fig. 1 . Homogeneous atomic chain 



The energy of a short bond is given by the formula 



while that of a long one by 



W f = -K f (V lU ) 2 + n f V lU 



W s = ^/c s (V 2 w) 2 + n s V 2 w, 



[I) 



(2) 



where Vi and V 2 are symbols of the central difference operators acting on displacements of 
the first and second neighbours, respectively. The elastic energy of the chain equals the sum 
over all the bonds. 

The state of equilibrium of a homogeneous chain is described by the set of equations 

Kf(Ui-i - 2Ui + U i+ i) + K s (Ui- 2 - 2ui + u i+2 ) = 0, (3) 
where the index i = 1,2,3,... . The general solution to equations (H) is given by the function 

Ui = Az l + Bz~ l + Ei + C (4) 
with arbitrary constants A, B, E, C, and 

Z=-P+y//P-l, (5) 

where 



P=- + l, a = ^. 
z, /Co 



(6) 



3 Atomic chain with "interface" 

Now, instead of a homogeneous chain we shall consider two semi-infinite homogeneous parts 
(of possibly different materials) connected by an 'interface' modelled by intermediate bonds 
of the third kind. The quantities related to the intermediate bonds will be labelled by the 
superscript °, while the superscripts ~ and + will refer to the left and right semi-infinite 
chains, respectively. 
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For simplicity we assume that the parameter a has the same value for all three materials, 

K -L= K J- = \ = a , (7) 

K~ K° K+ 

K g K, s K, s 

Moreover, we assume that 

a) the homogeneous chains are in equilibrium, 

b) there are no forces at infinity. 

The above conditions imply the following relation between the initial stretch forces: 

n f + 2n s = 0, (8) 

which means that the sum of forces transmitted by any section vanishes. 
The solution of such an interface problem is given by 

-^—z-i- x forj = -1,-2,... 
1 — z 

(9) 

A-^-z^ 1 forj = l,2,... , 
1 — z 

where z is defined by the formula (H). Displacements are calculated relative to the asymptotic 
value at — oo. The asymptotic value at +oo is given by 

A = s + £ -^±. (10) 

The symbols e_,e and e + denote deformations and are expressed by material parameters 
as follows: 

e^ = -T 1 (e + - o (n° s -n-)) , (11) 

e = 4(2/3 - T x - T 2 y\T x {n° B -nj)- T 2 (n+ - <)) , (12) 

e + = T 2 (e - i(n+ - <)) , (13) 
where we have introduced the auxiliary notations 

Tx = (1 - %1 + V 1 , (14) 

K° Z 

T 2 = + . (15) 

The resulting solution is represented by sharp-bend vertices of the line in Fig. 3. 

Let us notice that the number z is always negative. It implies a strongly oscillating 
character of the solution given by formula (|9|): the nearest neighbours are displaced in 
opposite directions. Thus every second atom is displaced in the same direction. Hence, as a 
result of the presence of the second neighbour interactions, we have obtained a polarization 
of the chain, expressed as a relative displacement of both substructures. 
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4 Composite atomic chain 



To construct the bi-continuum model we start from introduction of two fields of displace- 
ments: v and w, each of them defined on one of the sub-chains. 



b c I K f , n f K f , n f 

k r n 

\AAWAWvWWvWVVwAvvWWWWVVWAWAWWWVVV 



Fig. 2. Heterogeneous atomic chain 



We shall describe by the indices x and x the position of "black" and "white" atoms 
supporting the fields v and w, respectively. The description may be applied to two-component 
structures as well as to one-component but dimerized ones. Altogether, in the translationally 
invariant case we have four classes of bonds. Two bonds between the first neighbours are 
charcterised by parameters Kf, rif and Kf, fif, while two other bonds, between the second 
neighbours, are characterised by k s , n s and k s , n s . If there are no forces at infinity, then 
the sum of forces transmitted by any section vanishes, 



n s + nf + n s = 0, 
n s + n,f + n s = 0. 



(16) 



This fact implies, in turn, that initial stretches of bonds between the nearest neighbours 
equal each other: 

n f = n f . (17) 

The symbol a denotes now the distance between the closest atoms of the same colour (which 
corresponds to 2a in the previous sections), b and c denote the distances between atoms of 
different colours. These quantities are connected by the relation 



b + c = a. 



(18) 



The energy of the chain can be composed of 4 parts, each one connected with one kind of 
bonds, 

W = W S + W s + W f + W f , (19) 

where 
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w f = XX^M^ ~ w *) 2 + n A v * - w *))' 

Wf = 2 R A W * ~ V x-af + fif(Wx - V x - a )). (20) 

The last two parts of the energy contain terms combining values of different fields at different 
points. We shall express these terms with the help of fields v and w which interpolate the 
fields v and w to the whole chain. By definition, 

KfV x + K f V x _ a _ 

Vx = — -= , V x = V x (21) 

Kf + Kf 

and 

KfW x+a + KfW x _ 

w x = — — - — , w x = w s . (22) 

Kf + Kf 

Using these interpolating fields, we can express the energy of the first neighbour bonds by 
W f = ^2(-?; K f[( v x-Wx)-\ v— (w x+a -w x )] 2 + n f [(v x -v x )^ ^f— (% +a -%)]), (23) 

Kj + Aty Kf + Kf 

W f = y2(n K M w x-Vx) + — ^— (v x -v x „ a )f + n f [(w x -v x ) + Kf _ {v x -v x - a )]). (24) 

fif + fif fif + fcf 

All the differences are now expressed either by two fields at the same point, or by one field 
at different points. 

5 The bi- continuum 

We arrive at the corresponding bi-continuum expression for energy by appropriate trunca- 
tions of the Taylor series. As the simplest possibility let us consider the correspondence 

(1) Vf^af, 

(3) /— (25) 

The symbol / represents here an arbitrary function, the dash denotes the derivative with 
respect to x, and a - the lattice parameter. The set of rules fl25|) will be referred to as the 
correspondence of order 0. We shall consider also another set of rules, called correspondence 
of order 1, which takes into account the microstructure of the elementary cell of our composite 
chain. Namely, instead of rule (3) in (p5|) , we take 

(3a) v — > v H —v', 

Kf + Kf 

(36) w — ► w ^— w', (26) 

Kf + Kf 

where 

K = KfC — Kfb. (27) 
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If Rf = Kf and b = c, then both the correspondence rules are equivalent. Let us note that 
the transition to correspondence rules of order 1 does not augment the order of equations. 

Using the correspondence rule of order one obtains the following formulae for energies 
of the corresponding bonds: 

W s = J(^K s v' 2 + n s v')dx, (28) 
W- s = j \~R a w' 2 + n f w')dx, (29) 

Wf — — j [Kf{v — w H —aw') 2 + arif(v — w H — aw')]dx, (30) 

2(2 J Kf + Kf Kf + Kf 

Wf=^-f [R f (w - v H ^—av') 2 + an f (w —v-\ ^—av')]dx. (31) 

2d J Kf -\- Kf Kf -\- Kf 

In the last expression the equality ([[7]) has been used. The energy connected with the n.n. 
bonds may be rewritten in the form 

Wf + Wf = j-J [(Kf + Kf)(v - Wf + a 2 ^&_(KfW' 2 + KfV' 2 ) 

H J-^(Rfw' + Kfv')]dx + boundary term. (32) 

Kf + Kf 

The correspondence rule of order 1 introduce definite corrections to the above formulae. 
These corrections result in 

a) an additional coupling term of the form (v — w)(v' + w'), 

b) renormalisation of coefficients in other terms. 

In consequence, the energy density of the homogeneous bi-continuum takes the form 

Wd = ^ ev ' 2 + \s w ' 2 + k{v - w)(v' + w') + \h{v - w) 2 , (33) 

with the phenomenological coefficients. One can connect them with parameters of the dis- 
crete model by the following equalities: 

e = aK s + aKf' 



g = aK s + aKf\ 



Kf + Kf 
Kf + Kf 



Kf + Kf , 

h = i f , (34) 



and k is given by the formula (^7]). 

The equations of static equilibrium for a homogeneous medium have the form 

ev" — 2kw' — h(v — w) = 0, 

gw" + 2kv' + h{v -w) = 0, (35) 
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and, without forces at infinity, have the general solution 



Ae'^ + C, forx>0, 1 

Be^ x for £<0, 

Akte-^ + C, forx>0, 
w={ } (36) 

Bk 2 e^ x for a;<0. 

with the exponent coefficient expressed by material constants according to the formula 

^2 = He + 9)~ ^ 

eg 

The coefficients connecting amplitudes of fields v and w equal 

, h-C 2 e h-( 2 e 

h = - — , h = — — . 38 

n — 2Qk h + IQk 

In the case of exact inversion symmetry, when k s = K s , b = c, Kf = Rf and then k = 0, 

C 2 = — , fc! = A; 2 = -l. (39) 
In the next section we shall consider such a symmetric case. 

6 Comparison between bi-continuum and discrete so- 
lutions 



The solution of the bi-continuum interface problem may be constructed from functions (|3q), 
with appropriate values of constants A, B, C. The simplest way to match this solution to 
its discrete counterpart is to request that values of corresponding branches of functions u, v 
and w coincide at the atoms nearest to the interface. Two such conditions plus coincidence 
of asymptotic values will be sufficient to calculate A, B and C. In this way, making use of 
(f[2|), we have 

C = 7 M^I) (40) 



and, for the other constants 



A = B = -^£i-e Ca/4 . (41) 
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The bi-continuum fields v and w are plotted in Fig. 3. 




Fig. 3. Exemplary discrete and bi-continuum solutions; 
a = 0.05,K= 1.5, n= 1.5 



The upper branches correspond to the field w. The values of discrete function u correspond 
to the sharp-bends of the broken line. 
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Fig.4. Error in s.n.n. displacement vs. 1/oc 

To examine the accuracy of the bi-continuum model, the error in value of fields v and 
w at the distance 3a/ 4 from the interface, equivalent to the second nearest neighbours, is 
plotted vs. 1/a in Fig. 4. The error corresponding to a — 0.05 equals about 4 percent of 
the lattice constant. Note that, under the assumed correspondence rules, the error does not 
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change the sign. 

Another way to compare the atomic and bi-continuum models is to analyse the relation 
between the local dimerisation ranges or, equivalently, between the exponents in both models. 
Let us note that, since the symbol a denotes the distance between the second n.n., one should 
calculate the discrete exponent A from the equation 



- A a 



(42) 



which, according to (EJ), implies the dependence of A on a. On the other hand, for the 
bi-continuum model, by making use of equations ([ 



and (|34D, one obtains 
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a 



a 2 4 + a 



(43) 



The above equation implies 0<a(<A. 




Fig.5. Local dimerisation range, l/X vs. 1/^ 



The dependence of 1/A on is plotted in Fig. 5 as the lower curve to be compared with 
the line A = (. A substantial difference is present only at distances smaller than the lattice 
constant a. 
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To enable a more precise estimation of the error, the difference 1/A — l/( is plotted vs. 1/a 
in Fig. 6. 
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Fig. 6. Difference between dimerisation ranges: 
l/X- 1/£ vs. 1/a 



7 Final remarks 

We have presented a method of construction of the bi-continuum model based on a heuristic 
system - a linear atomic chain. One may, however, forget about the derivation, and consider 
the obtained equations as phenomenological ones. Generally, it is not necessary to relate 
the phenomenological coefficients to the parameters of an atomic model. The fairly good 
agreement between the bi-continuum and atomistic solutions can be treated as an argument 
in behalf of our model, at least in the sense that we have not made gross errors. 

The internal structure of the "atoms" in the chain may be complex - they need not be 
interpreted as material points. As an important example one can mention atomic planes 
of 3D crystalline structures, in particular the oxygen-copper planes, or even the octahedral 
layers in high temperature superconducing materials |§ . If our " atoms" are planes or layers 
||, one should replace, in the expression for energy, the scalar quantities by apropriate 
vectors and matrices. And the construction will run analogically, under the condition that 
one takes into account the interactions of the first and second n.n. between planes or layers. 
Similarly one can consider layered structures in the form of arrays of Josephson junctions. 
Such systems occur in a natural way in superconducting crystals, || or else they can be 
shaped artificially in a technological way. 

How important are the planar interfaces? There are many cases where they are important 
by themselves. It is expedient to notice that (what we call) the "planar" interfaces need not 
be strictly planar. It is sufficient that the curvatures of the corresponding surfaces should 
be not too high as compared with the characteristic exponents (both of dimension cm -1 , 
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inverse length). 

The interfaces can also affect the bulk properties of materials. It concerns the materials 
which contain a lot of interface area per unit volume. Also this property can be expressed 
by a parameter of dimension of inverse length which is particularly high for nanostructured 
materials. 
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